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Abstract 

Background: Pachysolen tannophilus is a non-conventional yeast, which can metabolize many of the carbon 
sources found in low cost feedstocks including glycerol and xylose. The xylose utilisation pathways have been 
extensively studied in this organism. However, the mechanism behind glycerol metabolism is poorly understood. 
Using the recently published genome sequence of P. tannophilus CBS4044, we searched for genes with functions in 
glycerol transport and metabolism by performing a BLAST search using the sequences of the relevant genes from 
Saccharomyces cerevisiae as queries. 

Results: Quantitative real-time PCR was performed to unveil the expression patterns of these genes during growth 
of P. tannophilus on glycerol and glucose as sole carbon sources. The genes predicted to be involved in glycerol 
transport in P. tannophilus were expressed in 5. cerevisiae to validate their function. The 5. cerevisiae strains 
transformed with heterologous genes showed improved growth and glycerol consumption rates with glycerol as 
the sole carbon source. 

Conclusions: P. tannophilus has characteristics relevant for a microbial cell factory to be applied in a biorefinery 
setting, i.e. its ability to utilise the carbon sources such as xylose and glycerol. However, the strain is not currently 
amenable to genetic modification and transformation. Heterologous expression of the glycerol transporters from 
P. tannophilus, which has a relatively high growth rate on glycerol, could be used as an approach for improving the 
efficiency of glycerol assimilation in other well characterized and applied cell factories such as 5. cerevisiae. 
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Background 

Pachysolen tannophilus is known for its ability to fer- 
ment D-xylose, one of the major components of hemi- 
cellulose plant residues, to ethanol [1]. However, 
P. tannophilus has also been shown to be capable of 
converting crude glycerol to ethanol under microaerobic 
conditions [2]. This ability is interesting since glycerol, a 
by-product of biodiesel production, has also been con- 
sidered as a potential alternative carbon source for in- 
dustrial bioprocesses due to the recent dramatic increase 
in production of biodiesel 

The whole genome of P. tannophilus CBS4044 has 
been sequenced [3], and with this, the possibility for 
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understanding and exploiting glycerol transport in this 
yeast has arisen. A number of studies have investigated 
glycerol transport, consumption and production by dif- 
ferent types of yeast [4-6], but the glycerol transport and 
metabolic pathways in P. tannophilus have not been 
studied so far. P. tannophilus has a relatively high 
growth rate on glycerol (0.3 h" 1 ) [2], but improvement of 
the strain for industrial applications through engineering 
the metabolism is hampered due to the lack of tools for 
genetic modification. However, with glycerol becoming 
an increasingly abundant substrate, it is highly relevant 
to characterize strains which are capable of its 
utilization, and consider the possibility of expressing 
genes from these strains in industrial cell factories for 
which genetic engineering strategies are already well 
defined. 
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Glycerol transport is the first barrier to glycerol 
utilization in a microbial cell In S. cerevisiae, glycerol en- 
ters the cell by two different mechanisms: a low affinity 
transport system (facilitated diffusion) and a high affinity 
proton symport system (active transport) [7]. When glu- 
cose is present, glycerol diffuses into the cells through a 
glycerol channel by facilitated transport, a process 
dependent on the FPS1 gene [8,9]. The physiological role 
of the facilitator Fpslp in S. cerevisiae has been de- 
scribed to be glycerol export rather than uptake during 
hypo-osmotic shock [9]. However, when only non- 
fermentable carbon sources (glycerol, acetate, ethanol) are 
present, an active uptake system driven by a proton 
motive force ensures the uptake of glycerol [7,10]. Two 
multi-membrane-spanning proteins encoded by GUP1 
and GUP2 (Glycerol Uptake Protein), were first identi- 
fied as being involved in active glycerol uptake in 
S. cerevisiae [10]. However, in later studies Guplp and 2p 
were proposed to have different roles than glycerol 
transport [11]. A screen for genes encoding membrane 
proteins involved in glycerol assimilation in S. cerevisiae 
identified a gene, STL1, involved in active glycerol up- 
take. Stllp is localized in the plasma membrane, is glu- 
cose repressed and inactivated by growth in glucose 
[7,12]. Importantly, it was concluded that the protein is 
a member of the sugar transporter family and acts as a 
glycerol proton symporter [12]. Stllp has also been 
shown to have a function in glycerol uptake in several 
other yeasts. It has been shown that in Candida albicans 
glycerol was actively transported into the cells by a pro- 
ton symporter encoded by the C. albicans STL1 [13]. It 
has also been reported that S. cerevisiae strains harbor- 
ing the STL1 gene from D. hansenii slightly improved 
their growth and doubling times on glycerol [14]. 

After glycerol is transported into the cells, two differ- 
ent routes have been identified in yeasts for further as- 
similation: a phosphorylation route and an oxidation 
route. Within the first pathway, glycerol is dissimilated 
by glycerol kinase encoded by GUT1 and then by 
glycerol-3-phosphate dehydrogenase encoded by GUT2, 
which is located at the outer surface of the mitochon- 
drial membrane [6]. Another pathway is catalyzed by 
glycerol dehydrogenase encoded by GCY1 followed by 
dihydroxyacetone kinase encoded by DAK1 and DAK2. 
In Schizosaccharomyces pombe, it has been reported that 
glycerol is utilized solely by this pathway [15,16]. In both 
pathways, glycerol is converted to dihydroxyacetone 
phosphate, which then enters glycolysis. The pathways 
involved in glycerol utilization differ in different yeasts. 
In some yeast strains, all four enzymes in both pathways 
are present, but only one pathway functions for glycerol 
dissimilation [4,6]. In S. cerevisiae, it has been shown 
that glycerol is degraded by the phosphorylation pathway 
and that mutants lacking one of the two genes are 



incapable of utilizing glycerol [17]. Although the fermen- 
tative pathway was also discovered to be present in 
S. cerevisiae [18], the function is unknown. The produc- 
tion of glycerol has two functions in S. cerevisiae: redox 
balance and protection against osmotic stress as a com- 
patible solute (osmolyte) [6,19]. Glycerol is commonly 
produced in the cytosol of yeasts from the glycolytic 
intermediate dihydroxyacetone phosphate. This com- 
pound is converted into glycerol in two steps that are 
catalyzed by glycerol-3- phosphate dehydrogenase (Gpd) 
and glycerol-3-phosphatase (Gpp) respectively. Each of 
these enzymatic activities has two isoenzymes, Gpdlp 
and Gpd2p, Gpplp and Gpp2p [20]. 

Although there have been several studies on glycerol 
metabolism in S. cerevisiae, many industrially used 
baker s yeast strains and laboratory strains grow poorly 
on glycerol. For example, the specific growth rate of the 
well-studied lab strain S. cerevisiae CBS 8066 under aer- 
obic conditions was found to be very low, 0.010 ± 
0.002 h" 1 , during shake flask cultivations in minimal 
medium containing glycerol as sole carbon source [21]. 
Strategies based on adaptive evolution, and overexpression 
or knock out of the genes involved in glycerol metabolism 
have been performed with the aim of increasing the host 
strains growth on glycerol and production of different 
value-added products have been pursued [21,22]. By adap- 
tive evolution, strains with growth rates of up to 0.2 h" 1 
have been achieved [21]. However, one possible rate- 
limiting step during glycerol metabolism, glycerol trans- 
port, was not addressed in those approaches [23]. 

In this study, a comparison of the genome sequences 
of S. cerevisiae and P. tannophilus has been performed 
to study the genes involved in glycerol transport and gly- 
cerol metabolism in these two organisms. Quantitative 
real-time PCR was performed to compare transcript 
levels during growth on glycerol. The genes involved in 
glycerol transport in P. tannophilus were heterologously 
expressed in a S. cerevisiae STL1 knockout strain to as- 
sess their function and to evaluate whether they could 
possibly contribute as transporter genes to improve the 
growth of S. cerevisiae on glycerol. 

Results 

Identification of genes in glycerol metabolism and 
analysis of glycerol transporters in P. tannophilus 

In order to find homologous genes related to glycerol 
transport and metabolism in P. tannophilus, the genome 
sequence was searched by BLAST with the sequences of 
genes known to have these functions in S. cerevisiae as 
queries. In P. tannophilus, two putative glycerol facilita- 
tor and two putative glycerol symporter genes were 
found and named PtFPSl, PtFPS2, PtSTLl and PtSTL2, 
respectively. In the glycerol consumption pathways, the 
genes with high similarity to S. cerevisiae were found 
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and named PtGUTl, PtGUT2, PtGCYl, PtGCY2, PtDAK 
(Table 1). In the glycerol production pathways, the genes 
found were named PtGPD and PtGPP. It was noticed 
that the genes PtFPS2 and PtGUTl were located close to 
each other in the P. tannophilus genome. With the aim 
of predicting the functions of the putative genes, bio- 
informatics tools Blastx and Blastp were used, and func- 
tional domain predictions were applied. The results are 
presented in Table 1. 

In order to investigate the phylogenetic relationship 
among putative glycerol transporters in P. tannophilus 
and their homologues in other yeast strains, an align- 
ment was performed with Ptfpslp, Ptfps2p, Ptstllp, 
Ptstl2p as well as with other predicted or published 
transporter proteins available from GenBank, Genolevures, 
SGD (Saccharomyces Genome Database) and CGD 
(Candida Genome Database). Unrooted phylogenetic trees 
are presented in Figure 1. 

For FPS1, the lengths of PtFpslp, PtFps2p, PpFpslp, 
YaFpslp and GlpF were relatively short compared to 
other Fpslp proteins. PtFpslp was shown to be 54% 
identical to PpFpslp, 45% to YlFpslp, 35% to ScFpslp 
and only 31% to GlpF. PtFps2p was shown to be 60% 
identical to PpFpslp, 49% to YlFpslp, 34% to GlpF and 
only 32% to ScFpslp. PtFpslp and PtFps2p were 
grouped together with YlFpslp and PpFpslp, while they 
were in a separate branch than ScFpslp. For STL1, the 
identity between PtStllp and ScStllp was low, 34%, 



compared to the 56% identity observed between PtStllp 
and YIPtStllp. 

Quantitative real time PCR 

qPCR was used to compare the levels of transcripts dur- 
ing growth on glycerol to those on glucose as the sole 
carbon source. It was found that the genes involved in 
glycerol transport and assimilation PtFPSl, PtFPS2, 
PtGUTl, PtGUT2, PtGCYl, PtDAK were up-regulated 
on glycerol to different degrees compared to on glucose, 
while the genes involved with glycerol production 
PtGPD and PtGPP were down- regulated on glycerol 
(Figure 2). Among the genes PtFPS2 and PtGUTl were 
most up-regulated, by a factor of 19.6 ± 2.3 and 17.6 ± 
1.3, respectively, on glycerol compared to glucose. The 
transcript levels of PtFPSl, another putative facilitator, 
were up-regulated by a factor of 5.4 ±0.9 on glycerol 
compared to glucose, while the putative glycerol 
symporter genes PtSTLl and PtSTL2 were expressed al- 
most at the same level on glycerol and glucose. PtGUT2, 
which is a putative mitochondrial glycerol phosphate 
ubiquinone oxidoreductase, was also up-regulated 4.2 ± 
0.04 fold on glycerol. 

Performance of S. cerevisiae stllA harboring different 
glycerol transporter genes 

The glycerol symporter gene ScSTLl, which is respon- 
sible for glycerol transport in the absence of glucose, 



Table 1 Genes potentially involved in glycerol metabolism in P. tannophilus 



Annotated function 



Gene name 
S. cerevisiae 



Gene name 
P. tannophilus 



Size (bp) 



Identity/Similarity 



Accession no. 



Transport 
Glycerol Facilitator 

Glycerol Symporter 

Consumption 
Glycerol kinase 
G3P dehydrogenase 
Glycerol dehydrogenase 

Dihydroxyacetone kinase 

Production 

G3P dehydrogenase 

G3-phosphatase 



ScFPSl 



ScSTLl 



ScGUTl 
ScGUT2 
ScGCY0031 

ScDAKl 
ScDAK2 

ScGPDl 
ScGPD2 
ScGPPl 
ScGPP2 



PtFPSl 
PtFSP2 
PtSTLl 
PtSTL2 

PtGUTl 
PtGUT2 
PtGCYl 
PtGCY2 
PtDAK 



PtGPD 



PtGPP 



1170 
972 
1728 
1905 

1848 

1998 

936 

933 

1767 



1314 



705 



35/49 
32/51 
34/51 
31/50 

53/67 
52/66 
54/73 
51/70 
42/58 



67/80 



33/53 



[GenBank:JQ481631] 
[GenBank:JQ481632] 
[GenBank:JQ481633] 
[GenBank:JQ481634] 

[GenBank:JQ481635] 
[GenBank:JQ481636] 
[GenBank:JQ481637] 
[GenBank:JQ481638] 
[GenBank:JQ481639] 



[GenBank:JQ481640] 



[GenBank:JQ481641] 



The potential genes were identified by performing a BLAST search with known 5. cerevisiae genes as queries. Sequences were deposited in Genbank and the 
accession no. is listed in the table. Sizes given for P. tannophilus genes. 
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Fps1 

KlFpslpAAQ01788.1 

KmFpslpAAS47031.1 

ZrFpslpAAQ16650.1 
ScFpslpCAA38096.1 
GlpFNP 4 18362.1 
— YlFpslpXP 504820.1 



100 



97 



100 



100 



36 



66 



10 



stn 



59 



PtFpslp JQ481631 

PtFps2 JQ481632 

PpFps 1 p XP 002494229. 1 



100 



100 



DhStlp DEHA2A12364p (Genolevures) 

ZrSTLlp ZYRO0E01188p (Genolevures) 
Scstllp Ydr536wp (SGD) 

KIStllp KLLA0F26246p (Genolevures) 

DhStlp CAG87598.1 
CaStllp orfl9.5753p (CGD) 



51 



100 



36 



76 



100 



-PpstllpXP 002493344.1 
-PtStl2p JQ481634 



9S 



-YlStllp YALI0C16522p (Genolevures) 
-PtStllp JQ481633 



Figure 1 Unrooted phylogenetic tree based on alignment of 
predicted protein sequences. The tree was bootstrapped 1000 
times. Accession numbers are presented after protein names. The 
origin of the protein sequences from databases other than GenBank 
are given in brackets after accession numbers. Ca, Candida albicans; 
Dh, Debaryomyces hansenii; Kl, Kluyveromyces lactis; Km, 
Kluyveromyces marxianus; Pp, Pichia pastoris; Pt, Pachysolen 
tannophilus; Sc, Sacchromyces cerevisiae; Yl, Yarrowia lipolytica; Zr, 
Zygosaccharomyces rouxii; GIpF, the aquaglyceroporin protein from. 
E coli. 



was knocked out. In order to validate the function of the 
glycerol transporters from P. tannophilus, all the pre- 
dicted transporter genes from P. tannophilus were heter- 
ologously expressed in S. cerevisiae stllA. Meanwhile, 
the glycerol facilitator Fpslp and glycerol symporter 
Stllp in S. cerevisiae were also expressed in S. cerevisiae 
stllA. To analyze the performance of the recombinant 
strains, growth in defined medium with 1% (v/v) glycerol 
as the sole carbon source was analyzed. The stllA strain 
expressing different transporter genes PtFPSl, PtFPS2, 
PtSTLl, PtSTL2, ScFPSl, ScSTLl and the S. cerevisiae 
CENPK 113-5D strain were tested under aerobic condi- 
tions in shake flasks. 

The stllA recombinant strain expressing the trans- 
porter gene PtFPS2 from P. tannophilus and ScSTLl 



from S. cerevisiae showed growth on glycerol after 96 - 
hours of cultivation, while the stllA recombinant strain 
harbouring PtFPSl, PtSTLl, PtSTL2 and ScFPSl did not 
grow (Figure 3). The strains expressing PtFPS2 and 
ScSTLl genes grew on glycerol with a |i max of 0.07 ± 
0.008 and 0.09 ± 0.005 and achieved a final OD of 8.1 
and 10.5 at 96 hours cultivation, respectively. In com- 
parison, the S .cerevisiae CENPK 113-5D strain grew 
very slowly with a |i max of only 0.02 ± 0.004 on glycerol. 

In terms of glycerol consumption, the recombinant 
stllA strains harboring the genes PtFPSl, PtSTLl, 
PtSTL2 and ScFPSl did not consume glycerol even after 
96 hours of cultivation. The S .cerevisiae CENPK 113-5D 
consumed only around 0.1 g glycerol. However, the re- 
combinant stllA strains harboring the genes PtFPS2 and 
ScSTLl assimilated glycerol much faster than the other 
strains and the final glycerol consumed for the two 
strains was 6.2 ±0.15 g/L and 9.8 ± 0.6 g/L glycerol, re- 
spectively (Figure 4). Although succinate was chosen as 
a buffer in all the experiments, it was not consumed by 
the strains. 

Discussion 

With regard to glycerol transport, two glycerol facilitator 
homologues PtFPSl, PtFPS2 and two glycerol symporter 
homologues PtSTLl, PtSTL2 were found in the genome 
of P. tannophilus. Based on qPCR results, all the gene 
homologues identified to be involved in glycerol trans- 
port and assimilation were induced on glycerol as car- 
bon source relative to the levels observed on glucose. 
The exceptions were the predicted glycerol symporters 
PtSTLl and PtSTL2, which were constitutively expressed 
on both glycerol and glucose as carbon sources. PtFPS2 
and PtGUTl (19.6 fold and 17.6 fold) were the two most 
strongly up-regulated genes, which might lead to the hy- 
pothesis that they are most relevant for glycerol assimi- 
lation in P. tannophilus. Moreover, the two genes are 
located closely together in the genome in agreement 
with a model where the two genes are co-regulated and 
contribute to the same overall process, i.e. glycerol as- 
similation. Next, we focused our efforts on understanding 
the function of the genes involved in glycerol transport as 
this function may be the rate limiting step for glycerol me- 
tabolism. In this paper the functions of the two types of 
glycerol transporters from P. tannophilus were addressed 
by transferring them individually to the well characterized 
yeast S. cerevisiae. 

The low affinity transporter Fpslp is a glycerol facilita- 
tor protein and belongs to the major intrinsic protein 
(MIP) family of channel proteins with six putative trans- 
membrane domains (TMDs). Fpslp is responsible for 
transporting water, small molecules like glycerol, urea, 
NH 3 , C0 2 or ions without consuming energy. The 
physiological role of the facilitator Fpslp in S. cerevisiae 
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PtFPSl PtFPS2 PtSTLl PtSTL2 PtGUTl PtGUT2 PtGcyl PtGcy2 PtDak PtGPP PtGPD PtTAF12 

Figure 2 Relative expression levels of glycerol metabolism related genes in P. tannophilus on glycerol compared to that on glucose. 

The estimation of relative expression levels was based on 2" AACT , where AAC T = (C T gene of interest - C T internal control) sample A - (C T gene of 
interest - C T internal control) sample B. C T represents the cycle number at which a sample reaches a predetermined threshold signal value for the 
specific target gene. All the experiments were performed in triplicate. 



was previously described to be glycerol export rather 
than uptake during hypo-osmotic shock and the Fpslp 
channel closed and retained the glycerol inside the cells 
in response to hyperosmostic shock [9]. An N-terminal 
domain 225 LYQNPQTPTVLP 236 and a C-terminal domain 
535 HESPVNWSLPVY 546 were found to have important 
roles in controlling Fpslp function. The N-terminus was 
required for closing of the channels and restricted trans- 
port through Fpslp. It was found that the rate of glycerol 



efflux was higher than that for uptake [9]. In agreement 
with this, we find in the present study that overexpression 
of ScFPSl does not suppress the glycerol transport defect 
of stllA strains as they do not grow on glycerol. Surpris- 
ingly, overexpression of PtFPSl in stllA strains allowed 
growth and assimilation of glycerol indicating that this 
homologous transporter is involved in glycerol uptake and 
recovers the stll A strain growth defect. The similarities of 
PtFpslp and PtFps2p to ScFpslp are 35% and 32% 
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Figure 3 Growth of recombinant S. cerevisiae strains expressing different glycerol transporter genes. The growth was tested in defined 
medium containing 1% (v/v) glycerol as the sole carbon source at 30°C in agitated flasks. Results represent the mean of at least duplicate 
experiments. 

V ) 



Liu et al. Microbial Cell Factories 2013, 12:27 
http://www.microbialcellfactories.eom/content/12/1/27 



Page 6 of 10 



12- 



10 



o 6 
o 

O 4 




PtFPS2_XY2 

ScSTL1_XY6 

CENPK 113-5D 

PtFPS1_XY1 

PtSTLl _XY3 

PtSTL2_XY4 

ScFPS1_XY5 



— i — 

20 



- 1 — 

40 



60 



- 1 — 

80 



100 



Time (Hours) 



Figure 4 Glycerol consumption of recombinant S. cerevisiae strains expressing different glycerol transporter genes. The growth was 
tested in defined medium containing 1% (v/v) glycerol as the sole carbon source at 30°C in agitated shake flasks. Results represent the means of 
at least duplicate experiments. 



respectively, but the homology is only restricted to the 
core of the protein with the six putative TMDs. ScFpslp 
(669 amino acid residues) is much longer than PtFpslp 
(389 residues) and PtFps2p (323 residues), the size of 
which is more similar to the E. coli GlpF (281 amino acid 
residues). The size difference is mainly due to the long 
hydrophilic N- and C- terminal domains of ScFpslp, which 
are absent in PtFpslp, PtFps2p and GlpF (Figure 5). By ex- 
pressing PtFPS2 in S. cerevisiae stllA strains, the facilita- 
tors from P. tannophilus increased the glycerol influx and 
glycerol consumption, presumably due to improved 
glycerol transport ability. Previously it was also reported 
that the glycerol transport was increased approximately 
2.5-fold in S. cerevisiae by introduction of the bacterial 
gene GlpF [9]. Therefore, we propose that in 
P. tannophilus the glycerol facilitator might function for 
glycerol influx rather than efflux. 

The high-affinity transporter Stllp was identified as 
the glycerol proton symporter in S. cerevisiae, which is a 
member of the sugar permease family of the major 



facilitator superfamily (MFS) [24]. It was demonstrated 
that the transcription of STL1 gene was significantly in- 
duced with glycerol as the sole carbon source, and the 
STL1 gene was subject to glucose repression based on 
microarray-based transcriptome analysis [25,26]. However, 
glycerol uptake by STL1 from C. albicans was not affected 
by the carbon source and salt stress [13]. In D. hansenii, it 
was shown that the active glycerol transport system was 
constitutively expressed and not subject to glucose repres- 
sion [27]. In agreement with these experiments, we showed 
that the glycerol symporter genes PtSTLl and PtSTL2 in 
P. tannophilus are constitutively expressed on glycerol and 
glucose based on qPCR expression analysis. However, the 
presence of the PtSTLl and PtSTL2 genes had no obvious 
effect on the physiology of S. cerevisiae, while the glycerol 
consumption and growth in strains that overexpressed 
ScSTLl was improved compared to stllA strains. The 
symporter PtStllp showed a low degree of sequence 
identities to ScStllp with 34%. However, the PtStllp from 
P. tannophilus exhibited 52% identity to the Stllp from 



240 250 560 ^ 570 

ptFps2/i-323 v:". L ".":v: "jTiIedeVi d- vp ^qgq^pTnkpyfgfgrff 

PpFps1p_XP_002494229. 1/1-326 QEPDFSA I EEELP VFTGS V S P L N K P Y F G F G K F F 

PtFps1p/1-389 HGGP ATG L EYPGDEEEKDEP I YRG LETPCNQ LD YG Y L I FF 

YaFps 1p_XP_504820. 1/1-322 LYTGDESPL NWP Y MG F D R I F 

GlpF_NP_4 18362/1-281 I G R H LPC D I C VflEE - 

KtFps 1p_A Y 338364/1 -564 K P K I j^fflQSHDSl S T Y ' Y Q G HESPLNLPLS A Y T D WF 

KmFps1p_A Y54 1009/1-572 < - K I flBfTCTflflBHIg = T Y I Y Q G H E S P L N L P L S A Y T DWF 

ZrFps 1p_A Y2995 13/1-658 VRPK i^^g, j-^, j A * ' ° ^m^m^U^Ki ° ' 

N-terminal domain C terminal domain 

Figure 5 Multiple sequence alignments with Fpslp from different yeast strains and E. coli. The alignments were done by CLUSTALW2, 
showing the N-terminal and C-terminal extension domains of Fpslp. 
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D. hansenii (DEHA2A12364p) and 56% to the Stllp from 
Yarrowia lipolytica (YALI0C16522p). It has been reported 
that Y. lipolytica can grow on glycerol with a |i max around 
0.3 h" 1 [28] and Pichia pastoris can grow on glycerol with 
a Umax 0-26 h" 1 [29]. P. tannophilus can grow on glycerol 
with (i max around 0.29 h" 1 [2], while S. cerevisiae grows 
relatively slowly on glycerol with |i max of 0.02 h" 1 (CEN.PK 
113-5D). Both the facilitator and symporter similarities 
among P. tannophilus, Y. lipolytica and P. pastoris were 
higher than that compared to S. cerevisiae. However, more 
protein sequence data with verified functions are needed 
to provide a definitive conclusion. 

Since glycerol transport might be the rate-limiting step 
for glycerol utilization, heterologous expression of the 
glycerol transporters from yeasts which have relatively 
high growth rates on glycerol could be used as the ap- 
proach for improving the efficiency of glycerol assimila- 
tion in S. cerevisiae. Improved glycerol transport has 
been demonstrated here with an increased glycerol con- 
sumption rate and growth rate under aerobic conditions 
with S. cerevisiae. 

Conclusions 

The current study demonstrates the function of the gly- 
cerol transporters from P. tannophilus. Our studies open 
new possibilities for further improvement of glycerol fer- 
mentation in industrial yeast strains with heterologous 

Table 2 Strains and plasmids used in this study 

Plasmid or strain Relevant characteristics 

Plasmids 
PUG6 
pSP-G1 
pXI-5 
pSH47 
pXI-5-PtFPS1 
pXI-5-PtFPS2 
pXI-5-PtSTL1 
pXI-5-PtSTL2 
pXI-5-ScFPS1 
pXI-5- ScSTU 
S.cerevisioe Strains 
CEN.PK 11 3-5D 
st/ 7 A 

PtFPS1_XY1 
PtFPS2_XY2 
PtSTL1_XY3 
PtSTL2_XY4 
ScFPS1_XY5 
ScSTU _XY6 



expression of glycerol transporters from the glycerol util- 
izing P. tannophilus. This study thus proposes a possible 
route to development of glycerol-based bioprocesses in S. 
cerevisiae. 

Methods 

Strains and plasmids 

The P. tannophilus strain used in this study was 
CBS4044. S. cerevisiae CEN.PK 113-5D was used for the 
construction of the STL1 knockout strain. The plasmid 
PUG6 [30] was utilized as the template for amplifying of 
the loxP-kanMX-loxP cassette. The integrative USER 
vector pXI-5 [31] was used in this study for constructing 
the expression vectors. Plasmid pSP-Gl [32] was used as 
a template for amplifying the bidirectional promoter 
TEF1IPGK1. Plasmid pSH47 [30] was used for excision 
of the kanMX marker gene. All plasmids were propa- 
gated in Escherichia coli strain DH5a. All the plasmids 
and S. cerevisiae strains used in this study are listed in 
Table 2. 

Medium and culture conditions 

For quantitative real-time PCR (qPCR) experiments, 
P. tannophilus was grown at 30°C in defined medium, 
containing (per liter) 0.67 g yeast nitrogen base w/o 
amino acids, 10 g succinic acid and 6 g NaOH as buffer 
system (initial pH around 5.6), 2% (v/v) glycerol or 2% 



Source or reference 



Plasmid with loxP-kanMX-loxP disruption cassette [23] 

2u-based URA3 plasmid with TEF1/PGK1 promoter [25] 

Integrative USER vector, with reusable URA3 Marker [24] 

Centromeric plasmid, URA3, PGAL1-Cre-TCYC1 [23] 

pXI-5 with TEF1/PGK1 promoter and PtFPSl gene This study 

pXI-5 with TEF1/PGK1 promoter and PtFPS2 gene This study 

pXI-5 with TEF1/PGK1 promoter and PtSTU gene This study 

pXI-5 with TEF1/PGK1 promoter and PtSTL2 gene This study 

pXI-5 with TEF1/PGK1 promoter and ScFPSl gene This study 

pXI-5 with TEF1/PGK1 promoter and ScSTU gene This study 

MATa MAL2-8 C SUC2 ura3-52 Peter Kotter 

CEN.PK 1 1 3-5D St/7 A (490,1 279 )doxP This study 

stllA with plasmid pXI-5-PtFPS1 This study 

sf/7A with plasmid pXI-5-PtFPS2 This study 

sf//A with plasmid pXI-5-PtSTL1 This study 

sf/7A with plasmid pXI-5-PtSTL2 This study 

sf/7A with plasmid pXI-5-ScFPS1 This study 

sf/7A with plasmid pXI-5- ScSTLI This study 
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(w/v) glucose as the carbon source. Cells were harvested 
at mid exponential growth phase at the same optical 
density by centrifugation at 5000 g, 4°C for 5mins, 
resuspended in 1 ml ice cold ddH 2 0, and the cell pellet 
was stored at -80°C S. cerevisiae stllA harbouring dif- 
ferent glycerol transporter genes were cultivated at 30°C 
for 96 hours in agitated flasks with defined medium, 
containing (per liter) 7.25 g synthetic complete (SC) 
powder, 10 g succinic acid and 6 g NaOH as buffer sys- 
tem, 1% (v/v) glycerol as the sole carbon source. 

Sequence retrieval and analysis 

The genome sequence from P. tannophilus CBS4044 
was obtained from previous work and is available in the 
EMBL database [3]. The sequences of genes involved in 
glycerol transport and metabolism in S. cerevisiae FPS1, 
STL1, GUT1, GUT2, GCY1, Dakl/2, GPD1/2, GPP1/2 
were used as queries in a BLAST search against the gen- 
ome sequence of P. tannophilus. The genes with high 
similarities and high identities were listed as potential 
orthologous genes. Gene sequences were registered in 
GenBank at NCBI. The multiple sequence alignments 
with the amino acid sequences of FPS1 and STL1 trans- 
porters from different yeast strains and E. coli were 
performed by using ClustalW2 free program at PDBe. 
Phylogenetic analyses were performed by using PAUP* 
4.0b 10 [33]. Unweighted parsimony analysis was performed. 
Trees were inferred using the heuristic search option with 
tree bisection-reconnection branch swapping and 1000 ran- 
dom sequence additions. Maxtrees were unlimited, branches 
of zero length were collapsed and all multiple parsimonious 
trees were saved. Clade stability was assessed in a bootstrap 
analysis with 1000 replicates. Trees were visualized in 
Treeview [34]. 

Quantitative real-time PCR 

Total RNA isolation and cDNA synthesis 

Total RNA was isolated from frozen cells using an 
RNeasy Mini Kit (QIAGEN, USA) according to the man- 
ufacturer s protocol. The quantity and quality of the iso- 
lated RNA were measured by NanoDrop ND-1000. The 
total RNA was treated with DNasel (Fermentas) prior to 
cDNA synthesis. Five ug total RNA were used to 
synthesize cDNA employing the RevertAid™ First Strand 
cDNA Synthesis kit (Fermentas) following the manufac- 
turers recommendations using oligo(dT) 18 primer. The 
cDNA was used as the template for quantitative real- 
time PCR for determining the transcript levels of the tar- 
get genes under different growth conditions. 

qPCR was performed on a Stratagene Mx3005P using 
the SYBR Green technology. The qPCR reaction mixture 
was prepared with 5 ul of 5 times dilution of cDNA as 
template, 10 ul SYBR Green master mix, 2 ul of 1 uM 



forward primer and 2 ul of 1 uM reverse primer and 
ddH 2 0 to 20 ul. The PCR program for qPCR was as fol- 
lows: 10 min of incubation at 95°C, followed by 40 cycles 
of 95°C for 30 s, 58°C for 30 s and 72°C for 30 s, and fi- 
nally the temperature was increased from 55°C to 94°C 
to check for unspecific products. The number of fluores- 
cence threshold cycles (Ct) was calculated with the set 
threshold value by using Mx3005P software. Results 
presented are mean values of three independent experi- 
ments. Suitable primer pairs for all the genes investi- 
gated were designed using the Primer Express software 
v3.0 (Applied Biosystems) software with the following 
parameters: product size 140 - 180 bp and melting 
temperature (Tm) 57-59°C. The primers used in this 
work are listed in Additional file 1. In addition, control 
reactions which include all components for qPCR, ex- 
cept for the reverse transcriptase, were performed to de- 
tect the genomic DNA contamination. The absence of 
products under these conditions meant the absence of 
enough genomic DNA for successful amplification. 

The relative expression levels were calculated approxi- 
mately based on 2" AACT , where AAC T = (C T gene of interest - 
C T internal control) sample A - (C T gene of interest - C T 
internal control) sample B. C T represents the cycle number 
at which a sample reaches a predetermined threshold signal 
value for the specific target gene. Relative expression data 
were normalized to the relative expression value of the 
housekeeping gene TAF12 in each sample, thus giving 
normalized relative expression for a target gene. 

S. cerevisiae STL1 knockout strain construction 

The construction of the STL1 knockout strain was done 
through the loxP-kanMX-loxP /Cre recombinase system 
[30]. The primers for amplification of loxP-kanMX-loxP 
fragment from pUG6 were: KanMX_fw_STLl (sequence 
5'- AAAGCTGAAAATAGAGGGTTGCTGGTCAATTT 
AGAAGGTTCCACAATTGCCTGAAGCTTCGTACGC 
TG -30 and KanMX_re_STLl (sequence 5'- AAGCGT 
TTGTTGATGCACGAACTTTCATTGATGCAATTTCT 
GGTGGGTATCTAGTGGATCTGATATCAC -30. The 
loxP-kanMX-loxP fragment with 50 bp homologous 
sequences in both ends to the gene STZJwas transformed 
into CENPK-113-5D [35] for obtaining the knockout 
strain by selection colonies on YPD supplemented with 
200 ug/ml Geneticin (G418). The resistant cells were 
transformed with pSH47 and selected on SC medium 
lacking uracil. Positive clones were purified and grown in 
liquid YPD overnight. The next day, the cells were washed 
once with distilled water and then incubated in YPG medium 
for 2 h in order to induce expression of the Cre recombinase. 
The culture medium was diluted and spread on YPD plates. 
The colonies were picked and restreaked on YPD plates 
supplemented with and without G418. Clones that were not 
able to grow on the G418 containing medium were the 
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knockout strains. The knockout strains were streaked on 
5-FOA (740 mg/L) plates to eliminate the pSH47 plasmid 
Yeast extract-peptone-dextrose (YPD) medium, synthetic 
complete (SC) medium, and SC lacking specific amino 
acids were prepared as described previously [36]. 

Constructs harboring different glycerol transporters 

In this study, six expression plasmids carrying the 
glycerol transporter genes from P. tannophilus and 
S.cerevisiae were constructed. The USER vector pXI-5 
was digested with AsiSI and then with the nicking endo- 
nuclease Nb.BsmI for making the AsiSI/Nb.BsmI USER 
cassettes [37]. The proofreading PfuX7 [38] was used for 
amplification of fragments PtFPSl, PtFPS2, PtSTLl, 
PtSTL2, ScFPSl, ScSTLl and bidirectional promoter 
TEF1IPGK1 with appropriate USER tails for insertion 
into the designated USER cassette AsiSI/Nb.BsmI. 
Primers used are listed in Additional file 1. PCR was 
performed with PfuX7 DNA polymerase according to 
manufacturers instructions. DNA mixtures were pre- 
pared from purified digested vector, glycerol transporter 
gene, and bidirectional promoter TEF1IPGK1, 5x 
Phusion HF buffer (Fermentas), and 1 U of USER en- 
zyme (Fermentas), adjusted to 10 \A by adding ddH 2 0. 
The mixture was incubated at 37°C for 20 min, followed 
by 25°C for 20 min. The 10 \A reaction mix was used dir- 
ectly to transform chemically competent E. coli DH5a 
cells. All the glycerol transporter genes were cloned 
under control of the TEF1 promoter. Constructs were 
named as: pXI-5-PtFPSl, pXI-5-PtFPS2, pXI-5-PtSTLl, 
pXI-5-PtSTL2, pXI-5-ScFPSl, pXI-5-ScSTLl (Table 2). 
The plasmids were verified by sequencing and digested 
with NotI (Fermentas) for 1 h at 37°C and then 
subjected to gel purification. Each linearized construct 
was homologously integrated into the S. cerevisiae stllA 
genome by transformation; six yeast strain lines were 
constructed as summarized in Table 2. Transformants 
were plated onto SC-Ura for selection. The S. cerevisiae 
transformants were restreaked on SC-Ura plates for sin- 
gle colonies, and the desired homologous integration 
was checked by colony PCR. Sequences and descriptions 
of primers are presented in Additional file 1. 

Analytical methods 

Growth was monitored by measuring OD 450nm with a 
Shimadzu UV mini- 1240 spectrophotometer (Shimadzu, 
Japan). Samples were taken periodically from the flasks and 
filtered through a 0.22 \iM syringe filter, and supernatants 
were preserved at -20°C for later HPLC analysis. Concen- 
trations of the substrate glycerol in supernatants were mea- 
sured by HPLC refractive index detector RID-10A using an 
Aminex87H column (Bio-Rad, USA). Separations were 
performed at 60°C, flow rate of 0.6 ml/min and 5 mM 
H 2 S0 4 as mobile phase. 



Additional file 



Additional file 1: Primers used for Semi-quantitative RT-PCR and 
USER Cloning. 
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